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Abstract 

Using the density matrix renormahzation group and a bosonization approach, 
we study a spin-1/2 antiferromagnetic Heisenberg chain with near-neighbor 
couphng Ji and frustrating second- neighbor couphng J2, particularly in the 
limit J2 >> Ji- This system exhibits both dimerization and incommensurate 
spiral spin correlations. We argue that this system is closely related to a 
doped, spin-gapped phase of the one-dimensional Kondo lattice. 

I. INTRODUCTION 

Relatively little is known about the behavior of frustrated antiferromagnetic quantum 
spin systems, in comparison with their unfrustrated counterparts. In general, frustration 
reduces antiferromagnetic correlations and the tendency towards Neel order. In the presence 
of sufflciently strong frustration, classical systems often develop non-collinear sublattice 
magnetizations. The classical antiferromagnetic Heisenberg triangular lattice, for example, 
has three different sublattices with magnetization directions in a plane at 120° angles. The 
behavior of the corresponding quantum system is still controversial. Quantum systems may, 
instead, dimerize in the presence of frustration, as exemplified by the exactly soluble one- 
dimensional Majumdar-Ghosh model |I| (see below). Yet another possibility is some sort 
of spin-liquid state, without dimerization or sublattice magnetization, examples being one- 
dimensional nearest neighbor systems and possibly the Kagome lattice. For a number of 
reasons, frustrated systems tend to be difflcult to study — for example, quantum Monte Carlo 
typically cannot be used because the frustration introduces a minus sign problem. 



1 



FIGURES 



Jl 



FIG. 1. The zigzag spin ladder. 



In this paper, we study a one-dimensional spin-1/2 Heisenberg system, with near- 
neighbor couplings Jl and frustrating next-neighbor couplings J2- This system can also 
be considered a two-chain lattice with diagonal, or "zigzag" couplings, as shown in Figure 
m, and we will refer to the system as the "zigzag" chain. The model is described by the 
Hamiltonian 

H = T^lJ^Si ■ + Jaf, ■ f,+i + JiS, ■ (f, + f,+i)], (1.1) 

i 

—* — > 

where Si and Tj represent the s = 1/2 operators on the 2 chains. We are particularly 
interested in the case where J2 > Ji- 

Besides a general interest in understanding frustrated quantum systems, there are sev- 
eral additional motivations for studying this system. As we discuss below, a bosonization 
treatment of a generalized version of the one- dimensional Kondo lattice shows a close re- 
lationship between the spin degrees of freedom of that system when doped, and the zigzag 
chain. The evidence we present below for a spin gap in the zigzag chain may also indicate a 
spin-gapped phase in a Kondo chain. Physically, zigzag arrangements of atoms are common 
and this model is appropriate for real systems such as SrCu02, studied in Ref. M. 



FIG. 2. The spin ladder. 

This model exhibits some similarities, but also marked differences from the standard 
two-chain spin-ladder shown in Figure ^. While inter-chain coupling is relevant in the 

ladder model, producing a gap which scales linearly (up to logarithmic corrections) for either 
sign of the coupling, it is marginal for the zigzag chain, producing a gap only for antiferro- 
magnetic sign, which scales exponentially with coupling, A oc exp(— constant ^1/^2)- The 
marginal nature of the inter-chain coupling in our model, is reminiscent of weakly coupled 
two-chain Hubbard systems |^,^. However, the renormalization group analysis of our model 
is both simplified by the absence of charge excitations and complicated by the presence of 
irrelevant intra-chain operators with coupling constants of 0(1). Unlike in the spin-ladder 
model, our numerical results for the spin-spin correlation functions indicate that dimeriza- 
tion is present for all Ji < Jic ~ 4 J2 and incommensurate spiral correlations are present for 
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all Ji < 2J2. 

From the numerical point of view, for a ID or quasi-lD system, this model is quite 
challenging. As mentioned above, quantum Monte Carlo methods are not useful. Because 
of the very long correlation lengths, exact diagonalization is not useful for the range of 
parameters in which we are interested. The model can be studied with the density matrix 
renormalization group |jlO| (DMRG), but it is still a difficult system for several reasons. The 
exceptionally long correlation lengths mean that very long systems must be studied — in some 
cases we have studied systems with thousands of sites. Furthermore, the system, particularly 
for large J2, has two properties which slow the convergence of DMRG with the number of 
states kept per block: 1) the correlation length is long, and 2) it is composed of nearly 
independent chains. The slower convergence of DMRG with nearly independent chains is 
known from the two-chain spin ladder case. For J2/J1 = 2, for example, it was necessary 
to keep about 400 states per block for an adequate calculation of the spin-spin correlation 
function. In contrast, similar accuracy can be obtained in the S = 1 Heisenberg chain 
keeping 50 to 80 states. The calculations presented here were only possible because of new 
developments in the DMRG algorithms, resulting in increased computational capabilities by 
nearly two orders of magnitude, which will be mentioned here but discussed in more depth 
elsewhere. 

Previous studies of this system include the exact diagonalization work of Tonegawa and 



Harada [11|, the recent work of Bursill, et. al. |12| using DMRG and a coupled cluster 



method, and the DMRG study of Chitra, et. al. [13], who considered a more general model 



which included a dimerization term. This work extends and improves upon the previous 
work in several ways. Particularly for the difficult J2 > Ji regime, we have been able to 
obtain improved numerical results, which are consistent with the previous calculations. We 
have presented the first field theory treatment of the zigzag chain in the J2 » Ji regime, 
as far as we know, and found general agreement with the numerical results. In addition, we 
point out a potentially important connection between the zigzag chain and the generalized 
Kondo lattice. 

This paper is organized as follows. In Section II, we discuss a bosonization and renormal- 
ization group approach for the zigzag chain. In Section III, we discuss the DMRG methods 
used in the numerical calculations. In Section IV, we present numerical results for a variety 
of properties. In Section V we discuss the relationship of the spin ladder and the zigzag 
chain to the half-filled and non-half-filled Kondo lattice, respectively. Section VI contains 
conclusions. 



II. BOSONIZATION AND RENORMALIZATION GROUP TREATMENT OF 

THE ZIG-ZAG SPIN CHAIN 

Quite different field theory treatments of the zigzag spin chain are appropriate depending 
on the ratio Ji/J2- For small values of J2/ Ji it is appropriate to treat the model as a single 
chain. Extensive discussion of this model may be found in Refs. |]T^-0. The conclusion is 



that the model is gapless for J2 < J2C ~ 0.241167 [|17| , [r8| . We briefiy review the conclusions 
here. The low energy effective theory is a free massless boson with SU(2) symmetry, or 
equivalently a level k=l Wess-Zumino-Witten (WZW) non-linear a-model. The uniform 
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and staggered components of the spin operators are represented in terms of the left and 
right-moving density operators pl^r and the SU{2) matrix field g of the WZW model as: 



Sj ~ (pl + pr) + constant (—l)''tr(a(y'). (2.1) 

As well as the WZW or free boson Hamiltonian, there is an additional marginal interac- 
tion controlled by J2: 

^int ^ (^2 - J2c)pL ■ PR- (2.2) 

This is marginally irrelevant for J2 < J2c', that is it renormalizes to zero at long length scales 
but produces logarithmic corrections to the simple scaling behavior of the free boson model. 
For J2 > J2C it is marginally relevant; that is it renormalizes to large values producing an 
exponentially small gap and inverse correlation length: 

A = vj^ oc e-constantj,/(j,-j,.)_ (2.3) 
This massive phase is spontaneously dimerized. The order parameter: 

d =< 52i„i ■ S2i> - < S2i ■ S2i+i >oc< tr^ > . (2.4) 
Since tr g has scaling dimension 1/2 it follows that: 

docVA (2.5) 

As J2/ J I is increased further the correlation length decreases and the field theory description 
becomes less accurate. At the point J2I Ji = -5 the exact groundstates become the simple 
dimer configurations of Figure as first pointed out by Majumdar and Ghosh |l[]. Here the 
correlation length is only one lattice spacing; a rigorous proof of a gap has been given in 
this case W9\. 



FIG. 3. The two dimer groundstates of the Majumdar-Ghosh model. 

As J2/J1 is increased further, we expect the correlation length to increase again. Even- 
tually, when J2 >> Ji, it is more appropriate to think of two spin chains with a weak zigzag 
inter-chain coupling, as shown in Figure |l|. At Ji = 0, we obtain two decoupled chains 
with vanishing gap and inverse correlation length. For small J1/J2 we expect the gap and 
inverse correlation length to be small so we may perturb around the point Ji = using 
a field theory approximation. Only the low energy degrees of freedom of each chain are 
relevant. We may represent these as in Eq. (|2.1| ), introducing two sets of fields gi, pi^i and 
PR^i where z = 1, 2 labels the two chains. The difference between the ladder and the zigzag 
chain becomes evident upon bosonizing the inter-chain coupling. 

In the ladder model we obtain a coupling of the alternating spin components: 
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(2.6) 



This has dimension 1, so we expect it to produce a gap proportional to J2 (up to logarithmic 
corrections associated with the marginal operators). To check explicitly that all degrees of 
freedom are gapped it is convenient to switch over to abelian bosonization notation. The 
SU(2) fields Qi can be written in terms of the free spin bosons, (psi as: 

Here is the field dual to dx'Yl{x'y\. [The subscript s for spin is redundant here but 

we keep it to distinguish these fields from the charge boson that arises in the discussion of 
the Kondo lattice in Sec. V.] Thus the alternating spin operators are: 

tr((7(7j) oc (sin V27T(f)si, — cos V27T(f)si, sin V2TT(j)si)- (2-8) 

Introducing the sum and difference spin bosons: 



(f)sl ± (f)s2 
^s,± = 7= , 



(2.9) 



V2 

Hinta becomes: 

^inta ^ - cos{V47T(f)s+) + cos(V47r0s-) + (1/2) cos(y47r0s-). (2.10) 

Thus we obtain decoupled Hamiltonians for the two fields 0^-1-. The interactions are expected 
to produce gaps for both bosons for either sign of J2, proportional to | J2I. pO|-p^ 



In the zigzag model the coupling of alternating spin components cancels, and we are left 
only with the marginal coupling of uniform components: 

^int« = 2^1 (PLi + Pm) ■ (Pl2 + PR2). (2.11) 

The term (/T^i ■ Pl2+Pri ■ Pri) does not renormalize to lowest order; we will assume it can be 
ignored, apart from a small velocity renormalization. The remaining Hamiltonian is Lorentz 
invariant. Note that the unperturbed spin Hamiltonian separates into 4 terms for left and 
right moving spin bosons of type 1 and 2. Thus we may interchange p^i with f)R2i so that 
the interaction becomes diagonal in the index i, labeling the two species of spin bosons. This 
interaction is then precisely two copies of the one which occurs for a single s=l/2 chain, 
Eq. (p^) . [This decoupling is spoiled by irrelevant operators.] Thus we conclude that, for 
ferromagnetic zigzag coupling, the interaction renormalizes logarithmically to zero and there 
is no gap. For antiferromagnetic coupling there is an exponentially small gap: 

A = vs/i oc e-constantJ2/Ji_ (2.12) 

We might again expect a broken discrete translational symmetry. Since the marginally 
relevant interaction couples left and right movers on the two chain, it is natural to expect 
that the order parameter is: 
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d=<S2r-{f2,-f2,+^)>, (2.13) 

as shown in Figure Note that, when we regard the model as a single chain with first and 
second nearest neighbor interactions, this is the same nearest neighbor dimer order parameter 
that was discussed above. Thus there is need be no other phase transition for ~ 0.25 Ji < 
J2 < oc. In two-chain field theory language, the order parameter is ti^agi) ■ tr (0^(72). This 
has dimension 1, so it should scale linearly with A 

d oc A oc e-^o^s*^^t-^^/'^\ (2.14) 



unlike in the other limit, discussed above where d oc v A. Note that, in this phase, the much 
stronger intra-chain correlations, < Tj ■ Tj+i > are translationally invariant; only the very 
weak inter-chain correlations break translational symmetry. 

:^ 
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FIG. 4. Dimerization in the zigzag spin chain. 

In addition to the long-range dimer order, there is also finite-range Neel order, with 
correlation length ^. ^ becomes very long when J2 is only slightly larger than ~ 0.25 and 
again when J2 » J\. Near J2C this long but finite range order is of standard Neel type: 

< Si ■ >oc (-l)^-^e-l'-^'l/^. (2.15) 

In the other limit, J2 » Ji, the spins on each chain exhibit finite range Neel order: 

< Si ■ Sj >oc< fi ■ fj >oc (-1)^-^6-1^-^1/^. (2.16) 

Note that from the single chain point of view, it is the second nearest neighbor spins which 
are Neel ordering. 

Further insight into the behavior of this model can be obtained from solving the classical 
problem. This gives a spiral order with pitch angle 6. 9 is the angle between neighboring 
spins in the single chain description. The classical energy per unit length is: 

E{e) = JicosO + J2Cos2e. (2.17) 

For J2 < O.25J1 the solution is the standard Neel state, but for J2 > Ji the pitch angle is 
given by: 

cos9 = -J1/AJ2. (2.18) 

This angle varies from tt at J2 — 0.25Ji to n/2 at J2/J1 — > 00. Of course, true Neel 
order is presumably impossible in a one-dimensional quantum antiferromagnet, even at zero 
temperature. However, we might expect the classical result to give a guide to possible 
finite-range order. Note that, in the case J2 >> Ji, the spins on different chains are almost 
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decoupled, | < S'j ■ Tj > | << 1, so classically we may regard these spins as orthogonal. This 
is consistent with the 7r/2 pitch angle predicted classically. We expect that the deviation 
of the pitch angle from 7r/2 defines a characteristic wave-length, which we expect to be 
proportional to i.e. 

0-7r/2ocl/^ (J2 » Ji). (2.19) 



For smaller values of J2I J\-, Chitra, et. al. |]T3[ showed (using DMRG) that the spiral 
correlations are only present for J2I J\ > 0.5, in agreement with the earlier work of Tonegawa 
and Harada []11|. Between the critical point and the Majumdar-Ghosh point, the pitch angle 



IS TT. 



III. THE DENSITY MATRIX RENORMALIZATION GROUP APPROACH 

The DMRG method is discussed in some detail in references iTDIl and lEBI. Here we 



discuss specifics of our calculations, as well as mentioning some recent improvements to 
DMRG which are implemented in our calculations. 

The numerical results were based on finite-system studies of systems with as many as 
600 sites, and on infinite-system results for systems up to 7000 sites. We obtain the total 
energies, bond energies, and the equal-time spin-spin correlation functions of the ground 
state. The accuracy of the calculations depend on the number of states m kept per block, 
as well as the value of J2/ Ji- In these calculations we used values of m up to 700. Large 
values of m were necessary only for large values of J2/ Ji- Near the Majumdar-Ghosh point, 
J2I J\ = 0.5, very small values of m sufficed; exactly at the Majumdar-Ghosh point, exact 
results for the ground state could be obtained with m = 2. Truncation errors, given by 
the sum of the density matrix eigenvalues of the discarded states, ranged from zero to 
O(10~^) for J2/J1 = 2.5. This discarded density matrix weight is directly correlated with 



the absolute error in the energy |TD|. We apply open boundary conditions to the lattice 
because the DMRG method is most accurate for a given amount of computational effort 
with these boundary conditions. 

We have been able to perform much more extensive calculations, both in terms of system 
size and number of states kept, than in previous DMRG studies. This is not because of 
more extensive computational resources-these calculations were performed almost entirely 
on a Digital AlphaStation 200 4/166 with 96 megabytes, rated at 135 SPECfp 92. The 
primary reason for the improved capabilities are some recent improvements to the DMRG 
algorithms. The main improvement involves a transformation of the wavefunction from the 
previous DMRG step to the current DMRG step, providing an excellent starting point for 
the sparse matrix diagonalization procedure. These improvements will be reported elsewhere 
In addition, the calculations were performed using a highly optimized C++ program, 
which, for example, translates all matrix operations into calls to a very efficient Basic Linear 
Algebra Subroutine (BLAS) library. This program keeps in memory only what is necessary 
at each step; everything else is written to disk. 
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IV. RESULTS 



The spin gap was calculated for J2/ Ji ranging from 0.4 to 2.0. The gap is defined as the 
difference in energy between the lowest 3^ = state and the lowest Sz = I state. For each 
value of J2/J1, this energy difference was calculated for a set of finite systems with open 
boundary conditions with different sizes L, with L typically ranging from 32 to 200. (We 
consider the system to be a L x 2 lattice, so the number of sites is actually 2L.) The gaps 
were extrapolated to L — 00 using a polynomial fit of the form 

AL = A + a2/L^ + a^/L' + ..., (4.1) 

and also including an ai/L term. The extrapolation was much more problematic than, for 
example, in the Heisenberg S = 1 chain, where the ai/L term should be omitted p5[. For the 



zigzag chain, the ai/L term appeared to be the dominant behavior except for exceptionally 
long systems. We dealt with this problem by using very long chains, so that the extrapolation 
itself was small. Our results are shown in Figure ^. 
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FIG. 5. The spin gap. 

We define the dimerization to be the absolute value of the difference of adjacent interchain 
bond strengths, given explicitly by Eq. ( p.l3| ). We calculated the dimerization using the 
infinite system DMRG method, extrapolating in the number of states kept per block m. 
The ground state is doubly degenerate, corresponding to a shift in the dimerization by one 
site. Ordinarily we expect to get a ground state with a fully broken dimerization symmetry, 
but it is possible numerical effects could produce an intermediate value for the dimerization. 
To prevent this, we add a very small dimerization field, proportional to Eq. (|2.13|) , to the 
Hamiltonian, with the coefficient taken as 10~^. A very long calculation was performed, 
starting with a relatively small m. Once the dimerization converged with system size for 
this value of m, m was increased, and allowed to converge again. Figure ^ shows the 
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most difficult case we studied, J2/J1 = 2.5. A maximum of m = 700 was used in this 
calculation, but substantially smaller values of m were necessary for most values of J2/ Ji- 
The plateau values are shown as a function of m in Figure |^. These can be fit very well with 
an exponential form, allowing the extrapolation to m ^ 00. (Near the Majumdar-Ghosh 
point, this procedure was not necessary: nearly exact results were easily obtained using 
small values of m.) The final results, corresponding to m ^ 00 and L 00, are shown in 
Figure ^ . 
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1000 2000 3000 4000 
/ 

FIG. 6. The dimerization during an infinite system DMRG calculation as a function of the 
DMRG step /. The numbers labeling each plateau are the number of states kept per block; once 
convergence was reached, the program was signaled to keep additional states. 
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FIG. 7. The dimerization plateau values from Figure 6. The solid line is a fit of the form 
y = 0.0405 + 0.03144 exp(-m/208). 
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FIG. 8. The dimerization as a function of J2I Ji- The solid line is a spline fit to the data, and 
the dotted line for J2/J1 > 2.5 is an exponential fit c? = 2.283 exp(—1.622J2/Ji). 

One interesting feature of these results is that the maximum dimerization does not occur 
at the Majumdar-Ghosh point, which one normally considers to be fully dimerized, with 
dimerization 3/4. A maximum of 0.7906135 occurs at approximately J2I J\ = 0.5781. (For 
the calculation of the maximum we reduced the value of the dimerization field to 10~^.) 
A dimerization greater than 3/4 is possible if one of the bonds is slightly ferromagnetic. 
Consequently, a maximum dimerization of unity would be possible if the two bonds were 
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independent. Since they are not independent, the maximum possible dimerization is re- 
duced. We can obtain the maximum possible dimerization by considering a model with the 
dimerization order parameter as its Hamiltonian 

H = Y.S2^■{f2,-f2^+l). (4.2) 

i 

The ground state of this system has the maximum possible dimerization. DMRG converges 
very rapidly for this Hamiltonian, yielding a dimerization of 0.8245165. Returning to the 
original zigzag Hamiltonian, we find that the weaker bond is ferromagnetic over a broad 
region: from 0.5 < J2/J1 < 2.5. The upper limit is somewhat uncertain, since we do not 
have reliable results for J2/J1 > 2.5. For large values of J2/J1, the dimerization can be fit 
very well by the function d = 2.283 exp(— 1.622 J2/ Ji), as predicted in Eq. ( p^.l4D , which is 
shown as the dotted line in the figure. 

The spin-spin correlation function C{i — j) = {Si ■ Sj) was calculated as a function of 
i — j using the finite system method, with i and j chosen symmetrically about the center of 
the system. The correlation function exhibits incommensurate behavior. This behavior is 
clearly seen in Figure ^ where we plot C{1) divided by the asymptotic form of the Lorentz 
invariant free massive boson propagator, to eliminate the dominant decay behavior of the 
function. 




50 100 150 
/ 

FIG. 9. The correlation function C{1) = {S ■ S), multiplied by /^/^e'/'», where the correlation 
length is ^ = 17.1. The system size was 200 x 2, and m = 350 states were kept per block. 

The correlation length ^ was chosen to make the maximum in the "beat" amplitudes as 
constant as possible. The equal-time free boson propagator is: 



dujdk 



(27r)2 a;7v2 + F + ^-2 



oc / dk 



VA;2 + e-2 



(4.3) 
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The fit was also performed without the factor of /^/^, with noticeably poorer results. Our 
results for the correlation length are subject to greater errors than, for example, the dimer- 
ization, both because a precise fitting function is not known and the errors in correlation 
functions are greater than in local measurements in DMRG calculations. 



The resulting values for ^ are shown in Figure |T0[ For comparison we show similar results 
for the ladder system shown in Fig. 2. The correlation length grows much more rapidly 
in the zigzag chain than in the ladder system. (These lengths correspond to considering 
the zigzag chain a L x 2 system. If we consider it a single chain, these lengths would be 
multiplied by two.) Figure |Tl] shows the same results as a semilog plot. It appears that for 
the accessible values of the correlation length does not yet grow exponentially with 

J2/J1. 
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FIG. 10. The correlation length as a function of J2I J\- 
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FIG. 11. A semilog plot of the correlation length as a function of J2I J\- 
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FIG. 12. The angle Q versus J2/J1. 
By fitting l^l'^^I^CiV) to a sinusoidal function with an arbitrary wavelength and phase, 



we were able to extract the incommensurate angle ^, which is shown in Figure |T2|. In Figure 
13, we plot ^ — 7r/2 versus 1/^, for which we expect linear behavior near the origin. The 
solid line corresponds to 



(20 



-1 



V/2 



1. 



(4.4) 
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Given the uncertainties in our procedures for extracting 9 and ^, our data points are rea- 
sonably consistent with this behavior. 




1^ 

FIG. 13. 61 - 7r/2 versus 1/^. 

V. KONDO LATTICE 

We consider a generahzed one- dimensional Kondo lattice Hamiltonian: 

H = Eh^(4c^+i + h.c.) + JS, ■ cf-c, + KS, ■ S,^^\. (5.1) 

i 

Here q^o- anihilates an electron of spin o at site % and Si is a spin-1/2 operator. Sums 
over spin indices are implicit. We include an explicit nearest neighbour spin coupling, K. 
This (and longer range terms) would be generated by conduction electron exchange (the 
RKKY mechanism). It could also arise from other exchange mechanisms in some cases as 
for instance in the organic chain compound CuPC(I). [^,^ We include it here because it 
allows for straightforward application of bosonization techniques. 

In fact the limit to which our method applies is J << t^K. In this limit, we are only 
concerned with the low energy degrees of freedom of the conduction electrons (q) and the 
spin chain. These can be represented in bosonized form, introducing spin (0si) and charge 
(0c) bosons to represent the conduction electrons and an additional spin boson (0^2) to 
represent the spin chain. This approach was attempted independently in Ref. but we 
disagree with their conclusions as explained below. The two spins bosons turn out to be very 
similar to the ones discussed for two spin chains above. They spin bosons may equivalently 
be represented by the matrix fields gi. For J = 0, we obtain simply 3 decoupled free boson 
Hamiltonians (plus various irrelevant interactions involving 0^2) • The velocities are vp = 2t 
for 0c and 0si and Vg = ttK/2 for 0^2- We note that the theory is not Lorentz invariant 
due to the difference of velocities. To bosonize the Kondo interaction we need the bosonic 
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representations for the conduction electron spin operator and the locahzed spin operators, 
Si. These are given by: 



2^^- ~ {PLI + pRi) + constant ■ [e^'^^^ti{agi)e''^'P^ + h.c] 

Sj ~ {pl2 + PR2) + constant (-l)-'tr (55(2). (5.2) 

Note that, away from half-fiUing, the continuum hmit Hamiltonian only contains a marginal 
coupling of the currents: 

H-mtu = JipLi + PRi) ■ {PL2 + PR2)- (5.3) 

At half- filling, the 2k f oscillation becomes commensurate with the alternating localized spin 
operator and an additional relevant (dimension 3/2) term occurs: 



^inta = ^tr(a5(i) ■ tT{ag2) cos(v27r0c), (5.4) 

where A oc J. 

The strong analogy with the two-chain spin system is now clear. The half-filled Kondo 
lattice resembles the standard spin ladder model, with a relevant inter-chain interaction. 
Away from half-filling the Kondo lattice is very closely analogous to the zigzag spin chain 
model. In this case there are no relevant charge interactions so we may expect a decou- 
pled massless charge sector. The field theory describing the spin sector is identical to that 
occuring in the zigzag spin chain, except for the difference of spin-wave velocities for the 
two chains. In fact, this situation could also easily be realised in the zigzag spin chain by 
having different coupling constants J2 and for the two chains. It is also clear that adding 
Hubbard interactions for the electrons does not change things very much, particularly away 
from half-filling. The charge excitations still decouple and remain gapless in that case. 
The scaling dimensions of various operators change with the Hubbard interaction strength, 
associated with rescaling the charge boson, in the standard way. 

At half-filling we expect H^^A^, which couples charge and spin modes, to produce a gap 
for all charge and spin excitations. The gaps should scale as for either sign of J. If 
we assume that the charge boson develops a gap and hence < cos(v^27r0c) >^ 0, then the 
interaction in the spin sector reduces to the same one that occurs for the spin ladder, Eq. 
( p.lOp . As argued above, this should gap both spin excitations. A different conclusion was 



reached in Ref. where it was claimed that does not appear in the interaction and 
hence remains massless. This incorrect conclusion was obtained because of a missing minus 
sign in the transformation from nonabelian to abelian bosons, Eq. ( p.7| ). The neccessity for 
the minus sign in the lower left matrix element can be seen by observing that the constraint 
detg = 1 is not obeyed and tT{ag) is not purely anti-hermitean, without it. (The "note 
added in proof" in Ref. pSf reflects a realization of this error. [p9|] ) 



Away from half-filling, where the charge boson is massless, we may analyze the spin 
sector much as for the zigzag chain. In particular, the Kondo interaction renormalizes to 
zero in the ferromagnetic case, leaving all spin excitations gapless. In the antiferromagnetic 
case, we expect an exponentially small spin gap: 
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Fujimoto and Kawakami p8|] assumed instead that J renormalized to some sort of strong 



coupling critical point corresponding to vanishing spin gap for 1 branch of spin excitations. 
This assumption seems rather unlikely from the point of view of the RG analysis, in light 
of the above comments, but was motivated by physical considerations. That is, it would 
seem that somehow the left-moving spin excitations from the localized spins interact with 
the right-moving spin excitations of the conduction electrons (and vice versa) to form a gap. 
While this seems reasonable at half-filling, it becomes difficult to understand away from 
half-filling. If we consider the strong Kondo coupling limit then localized spins form singlets 
with on-site conduction electrons. This clearly produces a gap at half-filling where there 
is one conduction electron for each localized spin. However, below half-filling there is an 
excess of localized spins which may produce gapless excitations. In fact, for the ordinary 
Kondo lattice model {K = 0) at strong coupling, J >> t, it has been shown by Sigrist et 
al. 0] that these left over spins form a (gapless) ferromagnetic groundstate. At weaker 



coupling Sigrist et al. found a non-ferromagnetic phase whose properties were not very well 



characterized. Fujimoto et al. ||2^ assumed that this phase corresponded to a single species 
of gapless spin excitations (as well as gapless charge excitations). 

We do not find this argument convincing and think that a spin-gap phase may occur at 
weak coupling. Evidence for this is provided by our analysis of the zigzag spin chain in the 
previous sections. Note that in that case also a transition to a gapless phase occurred for 
sufficiently large J, i.e. Ji > Jic ~ 4J2- 

A related phase with a spin gap but no charge gap has been found numerically in the t-J 
model away from half-filling, at J of order t. In this case, it is apparently not associated with 
any spontaneous discrete symmetry breaking, and may be thought of as a dimer fiuid state. 



3T| The absence of spontaneous symmetry breaking is a consequence of the gapless charge 
excitations. The dimer order parameter also contains a charge factor which has vanishing 
expectation value when the charge gap vanishes. 

There is actually a limit of the generalized Kondo lattice model, below half-filling, which 
is essentially equivalent to the t-J model: 

J » t,K » t'^/J. (5.6) 

The large J condition forces all conduction electrons to form singlets with localized spins. 
The unpaired localized spins can effectively hop around via the hopping term. Their pre- 
dominant interaction is the Heisenberg term, K. In addition they have various weak induced 
interactions of 0{t^ / J). These are the interactions responsible for ferromagnetism in the 
pure Kondo lattice model at strong coupling. In the spin-gap phase of the t-J model these 
interactions cannot change the behavior provided that they are small enough compared to 
the gap. This essentially constitutes a proof (given the numerical results on the t-J model) 
that the generalized Kondo lattice model has a spin-gap phase somewhere in its phase dia- 
gram away from half-filling. How large a region of parameter space is in the spin-gap phase 
and whether it includes the pure Kondo lattice model for some range of doping and J/t are 
open questions which we are investigating numerically. |^| 



It was argued independently by Zachar et al. that the Kondo lattice should have a 
spin gap away from half-filling using a different type of bosonization based on perturbing 
around a different limit of the model where the Kondo interaction is strongly anisotropic. 
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VI. CONCLUSIONS 



The zigzag spin chain is gapless for weak ferromagnetic interchain couphng, but has 
an exponentially small gap for small antiferromagnetic coupling. This phase has a weak 
spontaneous dimerization, or broken translational symmetry along with a finite-range in- 
commensurate magnetic order. 

Although we have presented results for a wide range of J2I J\-, our primary focus has 
been on the large J2I J\ region, where the system is best thought of as two weakly coupled 
chains. Most previous work has focused on smaller values of JijJx- Our results help to 
explain the behaviour of the quasi-one- dimensional antiferromagnet, SrCu02, studied in 
Ref. 0. This compound is believed to be well described by the zigzag spin chain (with 
very weakly coupled pairs of chains) with J2 ~ lOOOi^' and \ J2l J\\ in the range 10-1000. 
The susceptibility appears to go to a finite constant as T ^ 0, apart from a low T upturn 
attributed to impurities. This indicates the absence of a gap. This is to be expected 
from the results obtained here since the gap vanishes exponentially with Jij J\ and should 
be competely negligible for this range of couplings. Note that, if the interchain coupling 
had been of ladder type rather than zigzag type, this gap would have been much larger and 
perhaps observable in the available temperature range. We find that the gap is approximately 
.42 J' where J' is the inter-chain coupling in the spin ladder, for J' << J. Thus the gap 
might have been as large as 40K in the spin-ladder case and could have shown up in the 
susceptibility measurements which went down to 1.7K. 

The field theory description of the zigzag spin chain, in the limit, J 2 » Ji is closely 
related to the field theory description of the decoupled spin sector in the doped generalized 
Kondo lattice in the limit K,t » J. Our results on the zigzag spin chain suggest the 
existence of a spin-gapped phase in the doped Kondo lattice. 
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